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Ullemar's formula for the moment map, II 

Oh 

Abstract. We prove the complex analogue of Ullemar's formula for the Jacobian of 
the complex moment mapping. This formula was previously established in the real 



case. 



1. Introduction 

Consider the 'moment map' 

H : O -> (hq, fi u n 2 , ■ . .) 

where 

i 



Vk = — II Cd(Ad(, (i) 

and is a bounded domain in C If f2 is a simply-connected domain, we can uniformize 
it as the image Q = 0(D), where is a unique function which is holomorphic in the 
unit disk © and normalized by the following conditions 

0(0) = 0, 0'(O) > 0. (2) 



Then (1) takes the form 



M0) = 27 J J ^(*M*)l a «kAcfc. ( 3 ) 

D 

In general, when the function is not globally univalent in D and satisfies (2), we use 
the previous formula as the definition for the complex (or analytic) moments of [2]; 
then f2 is regarded CIS db Riemannian surface over C. 

This notion appears in several problems of complex analysis and its applications. 
In particular, the sequence (nk)k>i constitutes an infinite family of invariants of the 
Hele-Shaw problem [4] of the cell Q and can be used as a canonic coordinate system in 
the corresponding Laplacian growth model [3] (see also [6] for the functional analysis 
interpretation). 

In what follows, we consider the special case when the moment map fi is restricted 
to the set & n C C[z] of polynomials of degree n normalized by (2). It is easy to verify 
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that Hk — for all k > n (see also (13) below), so only the first n moments are of 
interest for P e & n . We consider the induced finite dimensional map 

fi: P = a lZ + a 2 z 2 + ... + a n z n -> (/j (P), ■ ■ ■ , A*„-i(P))- (4) 

Since /i (P) > 0, 

/i : S n -> R + x C n -\ 

Notice that 

dim M R + x U 1 - 1 = dim R & n = 2n-l, 

where & n is regarded as an open subset of a linear space. In [2] Gustafsson proved 
that the Frechet derivative d/i is non-singular at any P which is a locally univalent 
polynomial. 

On the other hand, notice that the subset @JJ which consists of the polynomials 
with real coefficients is an invariant set of \x in the sense that all /Xfc(P) are real (cf. 
(14) below). Hence, in a similar way the moment map induces the map //r : — > R n . 

In [7], C. Ullemar conjectured the following formula for the Jacobian of /xr: 

MP) ■= 9 ^-^\ P ) = 2-^ i a 1 ^P'(l)P'(-l)A(P-(,)), (5) 
0{ai, . . . ,a n ) 

where A stands for the principal Hurwitz determinant [1, §15.715], and Q* denotes the 
mirror conjugate image of polynomial Q, i.e. 

Q*(z) :=z m Q(l/z) =q m + q m ^ lZ + ... + q z m , m = degQ, (6) 



where Q(z) = Q(z) is the conjugate polynomial. 

The above expression for the Jacobian was recently proved in [5] as a consequence 
of the following identity 

J2 (P) = ^-l)™- 1 ^"^ K(P', P'*) ■ P'(-1)P'(1), (7) 
where TZ(-,-) denotes the resultant of the corresponding polynomials. 

In this paper we generalize formula (7) for polynomials with arbitrary complex 
coefficients. 

Theorem 1. The Jacobian of the moment map is expressed as follows 

J C (P) := g (/^-i,---,/xi^o^i,---,/Xn-i) = 2 a f-«+ 1 7^(P / ,P , *). (8) 
d(a n , . . . ,a 2 ,a 1 ,a 2 , ■ ■ ■ ,a n ) 

A well known theorem of Sylvester (see Section 2.2) allows us to compute the above 
resultant as the determinant of a matrix of size 2n— 2, whose entries are or a coefficient 
of either P' or P'*. In particular, the resultant is homogeneous in the coefficients of P' 
and P'* separately, with respective degree n — 1. 

On the other hand, geometrically, the hypersurface 

{(a,d): n(P',P'*) = 0} 
i.e. the critical set of the Jacobian, is the projection of the incidence variety 

n n— 1 

{(a,a,z): ^ka k z k ~ l = ^(n - k)a n ^ k z k } 

k=l k=0 
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that is to say, the set of (a, a) which appear above for some z. 

The following assertion is a direct consequence of the definition (15) of the resultant 
and formula (6) above, and it characterizes the set of critical points of dfi. 

Corollary 1. The moment map is degenerate at P if and only if the derivative 
P' has two roots oii and aj such that aidj = 1 (the case i — j is permitted). 

Note that for a locally univalent polynomial in the closed unit disk we have \ctj\ < 1 
for all the roots of its derivative. Hence, we obtain another proof of the above result 
due to Gustafsson [2]. 

Corollary 2. The moment map is locally injective on the set of all locally univa- 
lent polynomials in the closed unit disk. 

2. Preliminaries 
2.1. Complex moments. Using the Stokes formula, we obtain 

an an 

which implies 

M0) = 2< l + 1} / <t> k+ \z)$(z)dz = ±-. f <t> k {z)4>{z)<t>\z)dz, 

T T 

where T = <9B is the unit circle. Hence, using the identity z — 1/z which holds 
everywhere in T, we get 

M0) = 2ni( l +1) J 4> k+ \zW{l/z) $ = 2^/ fWfolMW dz - ( 10 ) 

T T 

Given a function which is analytic in a neighborhood of T, let us denote by A s (/) 
the sth Laurent coefficient of /, i.e. 

oo 

M = E A *(^ s > 

s=— oo 

hence 

l*<W = ^M<f +1 (z)4'0-/*)) = \-i{<P\z)<P\z)^l/z)) (11) 

Now, let P be an arbitrary polynomial in & n . Then P'(l/z) = P'*(z)z 1 ^ n and 
P(l/z) = P*(z)z~ n , which by virtue of (11) yields 

»k(P) = -^—\n(P k+1 P'*) = \ n -l(P'P k P*) (12) 

k + 1 

It follows from the first identity in (12) and P(0) = that 

fik(P) = 0, k > n (13) 
On the other hand, the second identity in (12) yields the so-called Richardson formula 

(14) 
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where the sum is taken over all possible sets of indices si, . . ., Sk > 1. It is assumed that 
dj = for j > n + 1. These formulae are easy to use for straightforward manipulations 
with the complex moments and it follows also that Hk{P) is a polynomial mapping. 

It is convenient to identify & n with the corresponding coefficient subset in R + x C n_1 
in a standard way: 

a ~ P := a\Z + a 2 z 2 + . . . + a n z a . 

Since, 

n 

Hq(P) = ^ s \ a s\ 2 > °) ^n-i{P) = na n x a n ^ 0, 

s=l 

the moment map (4) is well defined as an automorphism of & n into itself. 

2.2. Resultants. Here we review some basic facts about the resultant; see [8] for 
a detailed introduction. 

The resultant of two polynomials 

TO k 

A{z) = a m \\{z-aj), B(z) = b k ]J(z - (3 3 ) 



with respect to z is the polynomial 



U(A, B) = a k m bT J] (a t - ft 

i,j=l 



(15) 



The resultant vanishes iff A and B have a common root. It can be evaluated as the 
determinant of the Sylvester matrix, which is the following m + k by m + k matrix 

/a ai a m \ 



b 



a 1 
a 



bi 
bo 



ai 



bi 



a 
b k 



ai 



\ b h b k J 

in which the first k rows are the coefficients of A, the next m rows are the coefficients of 
B, and the elements not shown are all zero. The following are some useful elementary 
properties we will use below. 

U{A,B) = (-l) km n{B,A), 
U{A l A 2 , B) = n(A u B) n(A 2 , B), (16) 
K(z n ,A) = A n {0). 

Next, given a polynomial A(z) of degree n, we define its mirror conjugate image as 
A*(z) := z n A(l/z) = a n + a n _ x z + ... + a z n , 
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where A(z) = A(z) is the conjugate polynomial. We have for their roots: a* = (atj) 
and the corresponding resultant takes the following form 

f a ai a n \ 

Oq a l a n 



-1 



K(A,A*) = det 



On Q-0 

On On—1 ■ ■ ■ • • • 

On O n —i 



Oq 



\ 



Oq) 



(17) 



Remark 1. We wish to point out that the latter form, 1Z(A,A*), is irreducible as 
a polynomial of (a, a) over C. The proof is given in [5, Theorem 6]. 

3. Proof of the Theorem 

First, we evaluate the partial derivative of the moment map. Namely, we have for 
all k — 0, . . . , n — 1, j — 1, . . . , n, except for j — 1, k — 0, 

dfi k (P) 



X n ^(P'*P k ), 



da; 



(18) 



In fact, let j be an integer from {2, . . . , n}. Then by the first identity in (12) we 
have for 

dfi k (P) 1 



ddj k+1 



Xn(P 



fipk+1 

1 = \ n (P'*P k z j ) = \ n - j {P'*P k ) 



Similarly, using P* = a n + a n -\Z + . . . + a 2 z n 2 + a\Z n 1 and the second identity in (12) 
we obtain 



dfik(P) 



BP 

= \ n ^pip k °—) = \ n _ l (pip k z n -i) = \^{P'P k ). 



Finally, for j = 1 we have by the first identity in (12) 

%^ = r^7 A ^'*^ + P k+1 T~^ = *n- k (P'*P k ) + r^rMP^ 1 ). 

oai k+1 oai oai k + 1 

But Ai(P fc+1 ) = for k > 1, hence the desired assertion follows. 



We will make use the following notation 

df df df df df 



df df 



v <9a„' da n -\ ' da% da\ da% ' da n -i ' da n 
and by 

Qi-i = i a i 

we denote the coefficients of the derivative Q :— P' . Then 

V/i = (gn-i,---,5i,2go,gi,---,9n-i), 



(19) 
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and for all k — 1, . . . , n — 1 we have from (18) 

V// fc = (\ (Q*P k ), • • • , \n-2(Q*P k ), \n-l(Q*P k ), \l(QP k ), \n-l(Q*P k ))- (20) 

Let Y = V/io and for k > 1 write 

Y k := (a (q* A . . . , A n _ 2 (g*z fe ), A n _!(Q*^), Ai(g^), . . . , A n _!(g*^)). 

As a direct consequence of the above formula we conclude that 

Y k = 0, fc > n. 

Then it follows from (20) and 

P = z( ai + ... + a n z n - x ) 

that for all k > 1 

n-1 
j=k+l 

Thus, 

V/i A V/ii A • • • A V/i„_i = af Y A Y x A • • • A Y n _ u (21) 

where N = (n — l)n/2. 

On the other hand, for all k > 1 we have 

Y fc := (0, . . . , 0, g„_ b . . . , g fc , 0, . . . , 0, q , q u . . . , g fe _i), 

where the zeroes groups contain k and k — 1 items respectively. 

Now we treat the conjugate moments. We have fik(P) = f^k(P), whence 

V/Z fc = (VAi fc )*, 

where by X* we denote the mirror conjugate image of vector X = (xi,x 2 , ■ ■ ■ ,X2n-i), 
i.e. 

X* = {X2n-1, ■ ■ ■ , 

Repeating the above argument for the conjugate expressions yields 

V/Z„_i A V/i„_ 2 A • • • A V/ii = a? Y;.! A Y;_ 2 A • • • A YJ, (22) 

hence 

Vfi n -i A • • • A V/ii A V/i A V/ii A • • • A V;U n _i = 

= af' n Y;_! A • • • A Y* A Y A Y x A • • • A Y n _i. (23) 

We rewrite the latter identity in terms of determinants which gives the following 
expression for the Jacobian 
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where 





( 9o 






9n-l 














9i 


9o 




Qn-2 


q n -\ 












q n -2 


9n-3 




9o 9i 


q2 


93 • 




9n-l 




Y = 


q n -i 


q n -2 




9i 2g 


qi 


92 • 




9n-2 


9n-l 






9n-l 




q.2 9i 


qo 


9i • 




9n-3 


9n-2 










9n-l qn-2 








9o 


/7-i 

9i 




\ 






q n -i 










9o / 


and the elements not shown 


are 


all zero. 












Now, let Xj 


denote the jth column in Y. 


We have for 


j = 




,71 — 1 






X,= 


(o,. 


• -,0,5b, 9i, ■ 


-,Qn- 


i,0,... 


,0) 


T 




with j — 1 first zeroes, and for j 


= n + l,... 


2n- 


1: 












Xj = 


(o,. 


• • , 0, q n -i, ■ ■ 


•,9i,9o,0, ... 


,0) 


T 





with j — n first zeroes, and 

X n = (? n -i, • • • , 9i, 29o, 9i, • • • , 9n-i) T - 
One can readily verify that 

n-l 2n-l — 

Xn _J- q -^lX J + J2 — X j = (0, • • • , 0, 2g , 29i, ... , 29„-!) T , 

which yields for the determinant 
/ 9o 

9i 9o 9n-l 



det Y = 2 det 



\ 



9n-l 9ra-2 
9n-l 



V 



91 9o 9i 

92 9i 9o 

9n-l 9n-2 
9n-l 



92 

9i 



9n-2 9n-l 
9n-3 9n-2 

90 9l 

90 / 



The latter is the transposed Sylvester matrix of Q*(z) and zQ(z), hence by (16) 

det Y = 2 K(Q*, zQ) = 2{-l) n ^ K(zQ, Q*) = 
= 2 K(z, Q*) K(Q, Q*) = 2Q(0) K(Q, Q*). 

Thus, using our notation Q = P' we arrive at 

J C (P) =2af~ n+1 n(P',P'*), 

which completes the proof. 
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